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The Fermi -LAT has observed new 7 -ray flares from the blazar B0218+357 during July 2014. While no significant 
change in the 7 -ray spectrum has been previously observed through the flaring phase in late- 2012 , during this 
recent high activity the source displayed an exceptionally hard spectrum. The latter led to the detection of very 
high energy (VHE, E > 100 GeV) 7 -rays from B0218+357 by the MAGIC telescopes, establishing this source 
as the most distant TeV emitter known to date. In addition to the detection of VHE emission, this blazar is 
of particular interest since it is known to be a double-image gravitationally lensed system with a lens delay of 
11.46 ± 0.16 days measured in 7 -rays. We present the Fermi -LAT study of the July 2014 flares and discuss 
them in the context of previous measurements. 


I. Introduction 

B0218+357 is a blazar at z = 0.944 ± 0.002 [7j 
lensed by a galaxy at redshift 2 = 0.6847 [3]. The 
system appears as a double image in radio separated 
by 335 milli-arcseconds with a brighter western A and 
fainter eastern B images and an Einstein ring mm- 
The delay between the two images has been measured 
in radio as Af r = 10.5 ± 0.2 [2] and 10.1 ± 0.8 days 
[5|. Further timing analysis of the radio light curve 
from |6] using two independent methods resulted in 
two other possible delay values, A t r = 9.9+g or 

II. 8 ± 2.3 days 0. B0218+357 was detected for the 
first time in 7-rays by the Fermi Large Area Telescope 
[LAT; H| . Over the initial 4-years of LAT observations, 
B0218+357 displayed a 7-ray flux of F B>100MeV = 
(2.4 ± 0.9) x 10” 8 photons cm -2 s” 1 and a photon 
index of T = 2.28 ± 0.03 [3FGL J0221.0+3556; third 
Fermi- LAT Catalog, 3FGL. fl4l. 


2. Gamma-ray flares from B0218+357 

During the first six years of Fermi -LAT operations 
B0218+357 has displayed two major phases of en¬ 
hanced activity. The first was a long-lasting one that 
developed between 2012 June and 2013 March. It was 
characterized by a series of prominent 7-ray flares that 
persisted over four months. A detailed study of this 
period led to the measurements of the first robust 7- 
ray lensing delay. Applying an auto-correlation func¬ 


tion (ACF) to the evenly sampled light curve focused 
on the 2012-2013 enhanced 7-ray activity, a single 
dominant correlation peak was apparent. Its signif¬ 
icance was evaluated to be 9<r. The delay in the 7- 
ray data was inferred to be 11.46 ± 0.16 days (lcr). 
With the same dataset it was possible to estimate the 
flux ratios (A/B), which were found to be consistent 
with unity, and similarly, infer a magnification ratio 
of about unity J5]. Importantly, during this period no 
indication for spectral variability was found. The 7- 
ray photon index T = 2.30 ± 0.03 was consistent with 
the value derived before the start of the flaring; also, 
cf. the 3FGL value. 

B0218+357 underwent a second enhanced phase 
during 2014 July. The flares were between 2014 July 
11 and 27 (MJD 56849-56865) remarkably charac¬ 
terized by a distinctive hard spectrum |4]. Based 
on a preliminary analysis, on 2014 July 13 and 14 
(MJD 56851 and MJD 56852) the source was ob¬ 
served with respective daily averaged fluxes (E > 
100 MeV) of (6.5 ± 1.4) x 10” 7 photons cm -2 s” 1 
with T = 1.4 ± 0.1 and (6.7 ± 1.5) x 10”' photons 
cm” 2 s” 1 with T = 1.6 ± 0.1. Because the gamma- 
ray delay was previously measured, LAT Target of 
Opportunity observations were triggered to trace the 
expected delayed emission, about 11 days after. Note¬ 
worthy, emission at VHE was subsequently reported 
by the MAGIC telescopes at the time of the expected 
delayed 7-ray flare HI HU. In the following we present 
the preliminary analysis of the 2014 LAT dataset. 
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Figure 1: Fermi -LAT light curve computed in 6-hr bins at energies E > 300 MeV (lower panel) and the corresponding 
photon indices (upper panel). Vertical lines indicate the intervals of the flare (solid) and corresponding delayed emission 
(dashed). The two crosses indicate the arrival time of the two >50 GeV events (see Section [3.1| for more details). 


3. Hard spectral phase during 2014 

3.1. Fermi -LAT analysis 

Fermi -LAT data have been analyzed as described in 
[5| but restricted to the energy range 300 MeV - 300 
GeV to better characterize the hard-spectrum flaring 
behavior reported by [1] (ATel #6313). The source 
was modeled with a single power law with photon in¬ 
dex free to vary. In the analysis, 95% confidence level 
upper limits were computed when the test statistic 1 
(TS) [5] for the source was TS < 8 . A preliminary 
study of data between MJD 56849-56865 confirms 
the detection of B0218+357. During this interval the 
photon index is found to be on average exceptionally 
smaller (below 2) with respect to the value of 2.3 re¬ 
ported in previous LAT studies 0- Noteworthy, two 
high energy photons are positionally consistent with 
the source: a 95 GeV and a 54 GeV event, detected 
at ~MJD 56852, i.e. at the time of the initial flare. 

3.2. Light Curve 

Figure [l] presents the 6 -hr light curve (>300 MeV) 
and the corresponding photon index versus time vari¬ 
ations covering the time interval MJD 56849-56865. 


1 The TS corresponds roughly to the square of the signifi¬ 
cance assuming one degree of freedom. 


A first flaring phase, corresponding to flare of image 
A, is well defined by a fast flux increase accompanied 
by an exceptional hard spectrum and lasted 1.5 days. 
This phase is identified by six consecutive 6 -hr source 
detections, as denoted in Figure[l]by the two continu¬ 
ous vertical lines between MJD 56851.625-56853.125. 
Immediately after this short flare the source returned 
to quiescence, displaying only sporadic activity for the 
subsequent ~11 days. 


Then, at the time of the expected delayed emission, 
the source is observed again as denoted by the series 
of bins between the two dashed vertical lines, at MJD 
56863.125-56864.625. These detections happen ~11.5 
days later. The six consecutive 6 -hr intervals with TS 
> 8 are distinguished again by a small photon index 
(r < 2). This second flaring phase corresponds to the 
flare of image B, based on the hard-spectrum signa¬ 
ture. Formally, the ACF computed on the LAT light 
curve of 1-day binning between MJD 56849-56865 
shows a peak around 11 days (not shown). Though 
its statistical significance has not yet been estimated 
it is notable that the peak is broadly consistent with 
the previous 7 -ray delay of 11.46 ± 0.16 days. 


We compare the flux measured during the first 1.5- 
day flaring episode, beginning at MJD 56851.625, with 
the observation shifted by 11.5 days. We derive the 
resultant average flux ratio A/B ~ 4. 
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Figure 2: Scatter plot of the 6-hr binned fluxes and 
photon indices E > 300 MeV for the 1.5-day interval 
beginning at MJD 56851.625 (blue) and beginning at 
MJD 56863.125 (red), the same intervals reported in 
Figure [I] between solid and vertical lines, respectively. 

4. Summary 

Previous analysis of LAT data for the 2012 flares of 
the gravitationally lensed blazar B0218+357 found a 
7 -ray delay measurement of 11.46 ± 0.16 days and a 
flux ratio estimate between images (A/B) consistent 
with ~1. From a preliminary analysis of the 2014 July 
LAT data, we found: 

• Flares are characterized by peculiar hard spec¬ 
tra, with photon index T < 2 ; 

• Based on the hard spectral signature of the flar¬ 
ing emission, the delay is compatible with the 
previous y-ray delay measurement of 11.5 days; 

• The average flux ratio (A/B) over this flaring 
episode is ~ 4: while this factor is different with 
respect to the previous y-ray measurement, it 
is interestingly closer to the radio observations 

(~3 or4) m( 6 ]; 

• A 95 GeV and a 54 GeV photon were detected 
by the LAT during the first flare (A image) but 
were not observed during the second flare of im¬ 
age B. Nevertheless, the MAGIC detection en¬ 
sures that in the second phase the source emis¬ 
sion extends up to ^100-200 GeV. 
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